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ABSTRACT 

The WMAP team has recently highlighted the usefulness of combining the Ly-a 
forest constraints with those from the cosmic microwave background (CMB). This 
combination is particularly powerful as a probe of the primordial shape of the power 
spectrum. Converting between the Ly-a forest observations and the linear mass power 
spectrum requires a careful treatment of nuisance parameters and modeling with cos- 
mological simulations. We point out several issues which lead to an expansion of the 
errors, the two most important being the range of cosmological parameters explored in 
simulations and the treatment of the mean transmitted flux constraints. We employ a 
likelihood calculator for the current Ly-a data set based on an extensive 6-dimensional 
grid of simulations. We show that the current uncertainties in the mean transmission 
and the flux power spectrum define a degeneracy line in the amplitude-slope plane. 
The CMB degeneracy due to the primordial power spectrum shape follows a similar 
relation in this plane. This weakens the statistical significance of the primordial power 
spectrum shape constraints based on combined CMB-t-Ly-a forest analysis. Using the 
current data the simplest n = 1 scale invariant model with dn/d\nk = and no ten- 
sors has a Ax^ — 4 compared to the best fitting model in which these 3 parameters 
are free. Current data therefore do not require relaxing these parameters to improve 
the fit. 



1 INTRODUCTION 



Over the past several years we have seen the gradual emer- 
gence of a standard cosmological model. We seem to live 
in a universe dominated by dark energy and dark matter, 
with a small contribution from baryons. The universe is close 
to spatially flat and has adiabatic Gaussian initial condi- 
tions, consistent with inflation. This picture is conflrmed 
by a range of observations, especially those of the cosmic 
microwave background (CMB) anisotropies. The latest and 
most important confi rmation has been re cently reported by 
the WMAP sateUite dSpergel et al.ll2003h . 

Given that the classical cosmological parameters and 
the density content of the universe are now gradually be- 
ing settled the focus is shifting towards other cosmological 
tests that could reveal important information on the nature 
of the early universe and its creation. Among these is the 
amplitude and shape of the power spectrum of the primor- 
dial fluctuations. It can be constrained both by the CMB 
and other probes of large scale structure. Combining differ- 
ent data sets that cover a large range in scale is particularly 
powerful. Of the current cosmological probes, the Ly-a for- 
est - the absorption observed in quasar spectra by neutral 
hydrogen in the intergalactic medium (hereafter IGM) - has 
the potential to give the most precise information on small 
scales. 

For this reason, the WMAP team considered Ly-a for- 



est constraints in combination with CMB and galaxy clus- 
tering constraints to gain better leverage on cosmological 
parameters. The WMAP data alone are passably well fit by 
an inflationary CDM model with scale invariant spectrum 
n = 1. When combined with ACBAR, CBI and 2dFGRS the 
spectrum is pulled slightly towards n < 1 at A: = 0.05/Mpc 
and dn/dlnfc < 0. However, the evidence for running of the 
spectral index is only at the l.Scr level. When the data are 
combined with Ly-Q forest constraints the statistical signif- 
icance of both of these results is increased to 2-a. The pri- 
mary reason for this is that in the WMAP analysis of Ly-a 
forest the preferred slope at ~ 0.3 Mpc is n ~ 0.8. A run- 
ning spectral index at the level of dn/d\nk — —0.03 would 
have signiflcant implications for inflationary models. Here we 
will argue that the best estimate of the matter power spec- 
trum on Ly-a forest scales is higher in amplitude and slope 
than the WMAP analysis and that the uncertainty in the 
amplitude and slope of this power spectrum is larger than 
their estimate. Both of these changes weaken the case for 
a running spectral index and indeed make a scale-invariant 
n — 1 model appear to be a good fit to current CMB and 
Ly-a forest data. 

Ly-Q forest observations and constraints on cosmology 
have been explore d by several group s in the past after the 
pioneering work bv lCroft et alJ (|l99g)- On the o bservational 
side the two most recent analyses are those bv lCroft et alJ 
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J2002h . hereafter C02, and lMcDonald et alJ(l2000l) . hereafter 
MOO. The two groups obtain resufts for the flux power spec- 
trum, Ppik), in agreement with each other within the errors. 
In order to compare directly to the WMAP analysis we will 
only use C02 in this study. Both of these papers also explored 
theoretical implications. Addit ional recent theoretical ana l- 
yses have been performed bv [ G nedin fc Hamiltoiil i2002l) . 
hereafter GH, and Zaldarri aga et alTil20^E . hereafter ZSH. 
There are however significant discrepancies among them. For 
example, using a grid of PM simulations as opposed to di- 
rect inversion ZSIf show that the errors in C02 and GH 
may have been underestimated. Another apparent discrep- 
ancy is the different scaling of the linear amplitude of the 
power spectrum with the mean flux in C02 (equation 12) 
and GH (equation 8): the latter finds almost a factor of 5 
weaker dependence relative to the former. Yet another ap- 
parent discrepancy is the comparison between the C02 fit 
and t he WMAP reanalysis of the GH results iVerde et alJ 
l2003l) . While C02 conclude that a scale invariant n ~ 1 
model fits the data for Q.m ~ 0.3 and h ~ 0.7, the WMAP 
team uses the GH analysis of the same data (but the C02 
normalization uncertainty and expanded small-scale errors) 
to conclude that the Ly-a forest requires a shallow, n ~ 0.8 
slope on the scales probed by it (note that we are not re- 
ferring to the overall slope from the combined analysis of 
WMAP, just the Ly-a forest portion of it, see figure Ql. 

We will examine these differences in more detail below, 
here we highlight the most relevant theoretical issues: 

1) Coverage of cosmological models: the relation be- 
tween Ppik) and the 3-d linear theory mass power spec- 
trum, PL(k), is non-linear and model dependent, i.e., it has 
the completely general form Ppik) = f[PL(k)]. One must 
use cosmological simulations to connect the two. One can- 
not assume that their relation from one cosmological model 
applies to other models (at the level of accuracy needed for 
the present data) without verifying it with simulations. A 
grid of simulations must be built to cover the parameter 
space of interest. Some of the different scalings found in the 
literature may be simply due to expansion around a different 
fiducial model. 

2) Treatment of nuisance parameters, such as the 
temperature-density relation of the gas in the IGM and the 
mean transmitted flux, F{z) — (exp(— r)): as emphasized by 
C02, the inferred amplitude of the matter power spectrum 
is sensitive to the value of F. We will argue below that both 
the mean value and the error bar on F used by C02 were 
too low. Furthermore, the nuisance parameters couple non- 
linearly to the power spectrum and cannot be modeled as 
an overall calibration factor multiplying the inferred PL{k). 
Specifically, a change in F is not degenerate only with the 
overall amplitude of Phik) but also with its slope. This has 
been emphasized previously by ZSH and can be seen in figure 
16 of C02. However, both C02 and GH model the influence 
of F only as an overall calibration factor multiplying the 
inferred PL{k), ignoring its impact on the shape. 

3) Care with numerical details: the inference of the mat- 
ter power spectrum from Ppik) relies on numerical simula- 
tions. Previous investigations include some tests of the va- 
lidity of the approximate semi-analytic descriptions of the 
IGM used in these simulations and some investigations of 
the influence of numerical resolution and box size. How- 
ever, these tests do not demonstrate numerical convergence 



of the results at the level of precision that is of interest, given 
the quality of the flux power spectrum measurements them- 
selves. The numerical simulations used to interpret Ppik) in 
previous analyses were not tested for convergence of the final 
results with resolution and box-size. Similarly, GH and ZSH 
use collisionless particle-mesh simulations combined with 
approximate semi-analytic descriptions of the gas instead of 
fully hydrodynamic simulations, without quantitative tests 
of the accuracy of the approximations. 

The goal of this paper is to investigate the cosmological 
ramifications of a more reliable conversion between the mea- 
sured fiux power spectrum and the matter power spectrum. 
Of the three points mentioned above, all three are signifi- 
cant, but the combination of items one and two proves to 
be the most important for the current data sample. 



2 COMPARISON WITH PREVIOUS 
ANALYSES 

In the standard picture of the Ly-a forest the gas in the 
IGM is in ionization equilibrium. The rate of ionization by 
the UV background balances the rate of recombination of 
protons and electrons. The recombination rate depends on 
the temperature of the gas, which is a function of the gas 
density. The temperature-density relation can be parameter- 
ized by an amplitude. To, and a slope 7 — 1 = din T/d In p. 
The uncertainties in the intensity of the UV background, the 
mean baryon density, and other parameters that set the nor- 
malization of the relation between optical depth and density 
can be combined into one parameter: the mean transmitted 
flux, F{z). These parameters of the gas model. To, 7— 1, and 
F, must be marginalized over when computing constraints 
on cosmology. 

In observationally favored models, the Universe is ef- 
fectively Einstein-de Sitter at 2; > 2, so the cosmology in- 
formation relevant to the Ly-a forest is completely con- 
tained within Pt(fc) measured in velocity units. We pa- 
rameterize the linear power spectrum by its amplitude, 
A^(fcp) = kpPL(kp) , its effective slope, nofi(fcp) = 
dlnPi/dlnfclj. , and the effective running of the slope, 
dncff/dlnfc|^. , where kp = 0.03 s/km is the C02 pivot point. 

Rather than attempting to invert Ppik) to obtain the 
matter power spectrum we compare the theoretical Ppik) 
directly to the observed one. The advantage of this approach 
is that we can compute from the data without worrying 
about complicated window functions and covariances be- 
tween the data points (there will still be some covariances 
between the Pp{k) points, which are not given by C02, but 
simulations indicate that they will not be large on the rele- 
vant scales). 

For our present analysis we use a likelihood function 
module based on the C02 data and the library of simula- 
tions described in McDonald et al. (2003). These simula- 
tions cover the plausibly allowed range of F, To, 7 — 1, 
rieff, and dneff/dlnfc. Simulations with several box and grid 
sizes are used to guarantee convergence at 1% level, which 
is verified by detailed convergence studies on smaller box 
simulations. The grid is based on hydro-particle mesh simu- 
lations (Gnedin & Hui 1998; hereafter HPM), but these are 
explicitly c alibrated using fully hydrodynamic simulations 
iCen et alil994 . .2002,l . The simulation results are combined 
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in an interpolation code that produces Ppik) for any rela- 
tively smooth (CDM-like) input PL{k), F, To, and 7 - 1. 

We start by comp arinK C02, GH, and GH as interpreted 
bv lVerde et al.H2003h to each other. The original C02 result 
for the amplitude and slope of Pz,(fcp — 0.03 s/km) is plotted 
as the large filled square with error bars in figure a). Our 
fit to the points in Table 4 of C02 is plotted as the large 
open square with no error bars (some, but not all, of the 
discrepancy can be explained by our linearization of the fit 
using InP vs. Infc). C02 used a pure power law in their fit, 
and chose the pivot point to make the error bars on the 
slope and amplitude independent. We perform the fit using 
a ACDM transfer function shape, still using the power law 
free parameters, and plot the surprisingly different result as 
the open triangle in figure a). We suppress the error bars 
on the amplitude for this and subsequent fits to the C02 
or GH points because they are generally similar in size to 
the errors on the original C02 point. The difference in the 
form of the fit begins to explain the smaller n^ff implied 
by the WMAP analysis. Our identical fit to the points in 
Table 1 of GH produces the still small er n^ff represented 
by the small open square in figure nTal. lVerdeet alJ ll2003ll 
expanded the errors on the last 3 points to bring C02 and 
GH into agreement on these points within 1-cr. We show the 
results of fits using this prescription as the open pentagon 
and hexagon (for C02 and GH, respectively). Ultimately, the 
agreement between C02 and GH is not bad, but not perfect. 

Our version of the constraint is shown in figure 1(a) 
by the filled square for the best fitting point and the solid, 
68% contour for the error. This fit was performed assuming 
F = 0.705 ± 0.012 to match the C02 value (GH discuss a 
larger error bar but this only affects the error on the am- 
plitude - WMAP effectively used the C02 error). We will 
argue in Section 3 that this value and error bar are both 
too small. Our amplitude is 20% smaller than the GH am- 
plitude, while our slope is very similar. We disagree on the 
shape of the error contours - the GH errors on and Uef / 
are uncorrelated, while our contour shows significant correla- 
tion. The shape of our error contours is in a good agreement 
with those in ZSH, which also show a significant correla- 
tion between the slope and the amplitude. It also appears 
to agree with figure 16 of C02, although they did not bring 
this point out in their discussion of the Phik) errors. The 
errors one finds from a direct fit to the C02 o r GH points 
are to o small, but the expansion of the errors bv lVerde et alJ 
i2003h appears to have produced reasonable results. 

What is the reason for the disagreement between our 
central value, GH, and C02? We find that resolution and box 
size effects go in the right direction to explain the disagree- 
ment. GH used 20 Mpc/h simulations with 256"' particles, 
while the resolution used by McDona ld et al. (2003) is a fac- 
tor of two better than t his (see also iMcDonaldl i200ll) and 
iMeiksin fc Whitd j200lM . When we degrade our resolution 
to match GH we obtain the filled triangle and dashed con- 
tour in figure 1(a), with the increase in inferred slope and 
amplitude slightly over-shooting the GH result. C02 used 
simulations with another factor of two poorer mass resolu- 
tion. 

We also show the result of removing the hydrodynamic 
correction (the difference between hydrodynamic simula- 
tions an d HPM simulations ) that McDonald et al. (2003) 
(see also lMandelbaum et al.l (|20o3)) apply to the HPM sim- 
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Figure 1. Constraints on effective slope We// a-nd amplitude 
at kp = 0.03 s/km. (a) The large, solid square with error 
bars is the original C02 result for a power law fit, while the 
large empty square is our similar fit to the C02 Phik) points. 
Small empty symbols are our fits to the C02 and GH PL{k) 
results using a CDM shape of the transfer function (which dif- 
fers from a pure power law fit, which does not include the fact 
that the effective slope is rapidly changing even in the narrow 
range probed here) - triangle and square: fits to C02 and GH, 
respectively, using their error bars; pentagon and hexagon: sim- 
ilar with expanded errors on the last 3 points (slope errors are 
from the fit, amplitude calibration errors are all similar to the 
original C02 point). Small solid symbols and 68% contours are 
constraints from our direct fits to the C02 Pp(k) points — square 
and solid contour: our full procedure using C02 F constraints; tri- 
angle and dashed contour: similar with degraded simulation res- 
olution; hexagon and dotted contour: degraded resolution and no 
correction using fully hydrodynamic simulations, (b) The empty 
square and its associated errors show our estimate of the Ly-a 
constraint used by WMAP (we apply an 11% amplitude increase 
over the GH points). Solid symbols and contours are our fits to the 
C02 Ppik) points for different assumptions about F - hexagon 
and long dashed contour: F = 0.705 ± 0.012; triangle and short 
dashed contour: F = 0.705 ±0.027; pentagon and dotted contour: 
F = 0.742 ± 0.012; square and solid contour: F = 0.742 ± 0.027. 
The latter is our final result, for which we also show the 95% 
contour. 



ulations (filled hexagon and dotted contour). This has a 
comparable effect to the resolution and works in the opposite 
direction (note that this is only true for this specific choice 
of parameters). In the next section we find that the errors 
caused by these numerical details are insignificant compared 
to the current overall uncertainties, but they will become 
significant with future data sets. 
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3 MEAN ABSORPTION 

There are 3 nuisance parameters we include in our simu- 
lations: mean transmitted flux, F, mean temperature, To, 
and the slope of the temperature density relation, 7 — 1 
(we note that hydrodynamical simulations give results very 
close to the power law form). Among the nuisance parame- 
ters F has the largest effect, while the temperature-density 
relation leads to effects that can be neglected in the cur- 
rent discussion. There is considerable uncertainty in the 
value of F{z) [often transformed to an effective mean opti- 
cal depth, fsff = — ln(_F)]. There are two main approaches 
in the literature that attempt to deter mine it directly. The 
method pioneered bv lPress et al.Hl993ll . hereafter PRS, uses 
low resolution spectra and assumes an extrapolation of the 
quasar continuum from the red side of the Ly-a emission 
line. They give feff = A{1 + 2)"+^ with g = 2.46 ± 0.37 
and A = 0.0175 - 0.0056p ± 0.0002. Ignoring the errors one 
finds -fe// — 0.35 (here and in the following all the values are 
reported at 2 = 2.72, the central value for the C02 data). 
However, assuming the probability distributions for A and 
g are Gaussian one finds fe// = 0.29 ± 0.14. This value 
and error estimate appears to be inconsistent with their fig- 
ure 4, so it is not clear which to use. Note that C02 quote 
te// = 0.35 ± 0.018 and GH give f^ff = 0.35tg [;|4. Both of 
these, especially the former, have significantly smaller errors 
than one obtains from the original expression by PRS, but 
are closer (especially the latter) to a "by eye" estimate from 
figure 4 of PRS (note, however, that it is unclear how the 
band in PRS figure 4 was determined). 

A more recent analysis bv lBernardi et all (|20o3) uses 
a very similar method to PRS and applies it to a much 
larger SDSS sample, consisting of around 1000 quasars in 
the relevant redshift range. They find a rather similar result 
fe// ^ 0.35. While formal error bars are not quoted they ap- 
pear to be rather small, although the fact that the evolution 
they find is not very smooth complicates the error estimate, 
since one must be careful to match the redshift distribution 
to that of the C02 sample. 

The main systematic uncertainty in the PRS and 
iBernardi ct al. (2002J methods is the reliabil ity of the con- 
tinuum extrapolation: IBernardi et alJ ll2003ll assume that 
the continuum power-law found on the red side of the Ly-a 
emission line extends to the blue side, while PRS assume 
that the formula AX + BX^^^ can be used for this extrap- 
olation. In a recent analysis of low redshift HST quasars 
(for which the ab sorption is much le ss significant in the 
Ly-a forest range) iTelfer et al.l (|20o3) find that there is a 
clear break in the slope close to the Ly-a emission line, 
A ~ 1200 — I3OOA. The change in the slope is around one 
(see their figure 4 and table 1). This appears to produce a 
minimum 0.05 under estimate of F in the Ly -g forest region 
if one uses the PRS or lBernardi et al"] ll2003l) method. In this 
case instead of F = 0.70 one should use F = 0.75. 

A di fferent local c o ntinuu m fitting method has been ap- 
phed by iRauch et alj lll997t) and MOO to high resolution 
HIRES spectra, which have high signal to noise ratio. In 
this method one locally estimates the continuum using the 
regions where there is no apparent absorption. This is done 
by attempting to identify unabsorbed regions within the for- 
est and connecting them with a smooth curve. The main 
uncertainty in this method is the reliability of the local con- 



tinuum estimation, which would tend to overestimate F if 
there is some absorption everywhere. This method improves 
at low redshifts where the absorbers become rarer and the 
voids empty of neutral hydrogen can be identified more eas- 
ily. However, the number of QSOs analyzed was rather small, 
so the redshift evolution was only coarsely determined. In- 
terpolating MOO to z = 2.72 we find F = 0.742 ± 0.017, or 
fe// — 0.298±0.023, where the error is statistical only based 
on bootstrap resampling. This result is higher than the PRS 
value as given by C02 and GH. This was interpreted by C02 
as a sign of deficiency in the local continuum fitting method. 
On the other hand, the two values come into a close agree- 
ment if a 0.05 correction based on the Tolfor ct al. ( 200^) 
results is applied, arguing against the idea that the local 
continuum fitting method has a significant bias. 

Two other tests suggest that there is no large bias in the 
MOO F{z) estimate: First, to test the idea that the local con- 
tinuum method is biased by a lack of unabsorbed regions, 
MOO determined the mean value of the minimum absorp- 
tion in 10/i~^Mpc chunks of spectra from a hydrodynamic 
simulation. The offset, interpolated to 2: = 2.72, was only 
AF — 0.01, which should represent an upper limit on the 
bias. Second, we recently analyzed a sample of 15000 QSO 
spectra from the Sloan Digital Sky Survey to determine self- 
consistently the evolution of F{z) with redshift, the com- 
posite QSO spectrum and the principal components of the 
deviations from the mean composite (McDonald et al. 2003, 
in preparation). Since we do not assume anything about 
the shape of the continuum we can only determine F{z) as 
a function of redshift up to an overall constant, which is 
degenerate with the QSO shape. We find that the redshift 
evolution of MOO agrees well with this analysis within the 
errors, suggesting again that the systematic errors in the 
local continuum method are not large. If there was a bias 
from the local continuum fitting, one would expect the error 
to increase with redshift due to the increase in the mean 
absorption. That continuum extrapolation method of PRS 
may le ad to an underes tima te of F has also been pointed 
out bv lKim et all J200ll) and lMeiksin et al.l ||200J) . 

We should mention that C02 also present the fiux fill- 
ing factor statistic, which supports the lower value of F. 
This is a statistic based on the unsmoothed data and is sen- 
sitive to the smallest structure in the forest. As such it is 
probably sensitive to the resolution and small scale physics 
details in the simulations, as well as details of the noise 
distribution. Another promising method to determine the 
mean flux indirectly may be offered by the bispectrum anal- 
ysis, which concentrates on the large scales where both sim- 
ulation resolution and no ise properties are less important 
jMandelbaum et al.ll20o3) . 

Given the current situation it seems clear that using 
F — 0.705 ± 0.012 underestimates the error and probably 
the amplitude and that a more conservative treatment is 
needed. Here we wiU adopt the value F = 0.742 ± 0.027. 
We have slightly expanded the errors relative to MOO to ac- 
count for the uncertainties in the local continuum method, 
the redshift evolution, the possibility of imperfect removal 
of metal absorption by MOO, and the effect of damped Ly-Q 
systems. While we believe this value represents a more real- 
istic estimate, we will also show below the results assuming 
F — 0.705 but with expanded errors to more realistically 
account for the current uncertainties. 
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Figure^b) shows the change in the amplitude and slope 
of the inferred linear power at the pivot point if one uses dif- 
ferent values for F or for its error. We show both the best 
fitted value and the corresponding 68% error contours. The 
main effect comes from changing F from 0.705 to 0.742, 
even if we keep the errors unchanged at 0.012. We see that 
this changes the best fit amplitude by almost a factor of 2 
and the slope by 0.2. The effect on the amplitude is much 
larger than the GH estimate, but closer to the C02 estimate 
(we remind the reader here that the error estimates for the 
mean flux were smaller in C02 than in GH, making their 
final estimates comparable to each other and significantly 
smaller than our estimate). As discussed above, a change in 
F does not correspond to a simple rescaling of the overall 
amplitude of the matter power spectrum. One can see how 
the error contours really explode in the direction of large 

and large n^ff, probably a consequence of the fluctu- 
ations becoming more nonlinear. We note in this context 
that there are competing effects when one increases the am- 
plitude of fluctuations: while density fluctuations increase so 
do velocities and the corresponding suppression of power on 
small scales ("flngers of god" effect). As a result, the flux 
power spectrum can go up (on large scales) or down (on 
small scales), with the zero crossing not too different from 
the pivot point of the current data. The actual dependence 
of the flux power amplitude on the linear amplitude is thus 
rather sensitive to the range of scales covered by the data 
and the model around which one studies the dependence. 
This effect probably explains some of the differences found 
in the scalings between the linear amplitude and the mean 
flux. It highlights the need to use a grid of simulations that 
covers the parameter space of interest. 

In addition to the change in the mean value we also 
investigated the effects of increase in the error estimate. If 
we use F = 0.705 ± 0.027 the 68% contour again expands 
signiflcantly (figure [Tt). If we use F = 0.742 ± 0.027, our 
current best estimate, the error contours expand relative 
to -F = 0.742 ± 0.012, although the effect is less dramatic. 
The change in F moves the points and error contours along 
a degeneracy line diagonal in the A^-Ucft plane. Note that 
one cannot parameterize the error contours with a Gaussian 
form. 



4 COMBINED CMB-LY-a FOREST ANALYSIS 

To Eissess the effects of our improved treatment we present 
a joint WMAP+CBI+ACBAR ("WMAPcxt") - j-Ly-g for- 
est an alysis (I Verde et al. 2003; Mason ct al. 20(3 IKuo et aO 
l2002l). which can be compared to the WMAP results 
JSoergel et ai]l2003t |Peirise£aO Eooi) . While the WMAP 
team also included the 2dFGRS galaxy clustering results 
we have chosen not to do so in the present analysis. The 
main reason is that the 2dFGRS results do not significantly 
strengthen the significance of the results on the primordial 
power spectrum. While a joint analysis with the Ly-a for- 
est gives a 2-a indication of n < 1 and dn/dlnk < 0, the 
result with 2dF is actually less signiflcant than the result 
from the CMB alone and both n — 1 and dn/dlnk — 
are statistically acceptable at the 1.5-a level. In addition, 
the galaxy power spectrum is measured at low redshift, so 
one must include growth factor effects due to the possible 



nonstandard dark energy component. Finally, there remain 
concerns about the accuracy of the linear bias and redshift 
distortion assumptions on scales used in the 2dF analysis, 
as well as the remaining effects of luminosity bias. 

Our approach to p arameter estimation is a s tandard 
one. We run CMBFAST dSeliak fc ZaldarriaeallOQ^) to gen- 
erate both the CMB power spectrum and the matter power 
spectrum at z = 2.72. We use the former to obtain the CMB 
X^. We interpolate the latter onto the grid of simulations to 
obtain the for Ly-Q forest, marginalizing over all the nui- 
sance parameters (F, Tp, and 7 — 1). We generat e a Monte 
Carlo Markov chain JChristensen fc Meveill200lh using the 
combined x^- We obtain a sequence of chain elements that 
sample the probability distribution in the space of cosmo- 
logical parameters. We ran several chains both with and 
without tensors. We verifled that we obta in results similar 
to those reported by ISoergel et al.l diooi) using the CMB 
data only. 

It is instructive to look at the constraints in the A'^-UeS 
plane at the k — 0.03 s/km pivot discussed before. These 
are shown in figure (Altogether with our Ly-a contours. The 
CMB data define a degeneracy line in and ricfi, which 
is essentially determined by the curvature in the primor- 
dial power spectrum: low A^ and UcB indicates a very nega- 
tive curvature, while for high values of A^ and one has 
dn/d\nk > 0. Note that this degeneracy is very similar to 
the F degeneracy in the Ly-a forest. The filled triangles in 
figure 121 show the models produced by the CMB only chain 
with dn/dlnk > 0. One can see that while these models are 
between the 68% and 95% contours for the CMB, they are 
inside 68% for the Ly-a forest. In contrast, these models are 
2-(j away from the Ly-a constraint as used by the WMAP 
team (with or without tensors). The Ly-a constraint used by 
the WMAP team reinforces the slight preference for negative 
dn/dln k that appears in the CMB data alone. Our reanaly- 
sis of the Lya forest, by contrast, pulls away from the direc- 
tion of the best-fit CMB-only model towards dn/dlnk ~ 0. 

We ran Markov Chains with a combined 
WMAPext-|-Ly-a analysis with and without tensors 
(we fix the tensor slope to T/S = —Sut)- We only focus on 
the parameters that may be affected by our reanalysis (the 
primordial power spectrum and tensors) and we emphasize 
that the other constraints as presented by WMAP remain 
unchanged. The MC chain with tensors cannot be used 
to infer whether the data require running of the slope. 
This is because we find the running of the slope to be 
strongly correlated with the tensors, in such a way that a 
larger tensor contribution requires more negative running. 
This is not surprising, since negative running reduces 
power on large scales (compared to no running), which 
can then be filled in by tensors. It however suggests that 
any one-dimensional distribution of dn/dlnk will depend 
on the assumed prior for T/S and that positive dn/dlnk 
may appear unlikely only because of the bound T/S > 0. 
In the limit of a perfect degeneracy between the two 
parameters and assuming a uniform prior on T/S > the 
probability for dn/dlnk > simply decreases with the 
assumed upper cutoff on T/S. This means that adding two 
or more additional parameters when they are not required 
by the data can lead to misleading results, especially if only 
one-dimensional projections are used in interpretation. One 
way to address whether this is an issue is to compare the 
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Figure 2. Our Ly-a constraints (solid square and solid con- 
tours), CMB (WMAP+ACBAR+CBI) constraints (dashed con- 
tours) and the Ly-a constraint used by WMAP (empty square 
with errors). The CMB contours follow dn/dlnk, with more neg- 
ative values leading to low amplitude and low slope: very roughly, 
^cft ~ —2.45 corresponds to dn/dlnk = 0, n^g ~ —2.7 to 
dn/dlnk = -0.03 (final WMAP value) and n^g ~ -2.85 to 
dn/d\nk = -0.05 (WMAPext value). We warn that due to a 
small number of Markov Chain elements (around 25000) the con- 
tours shown here may be an underestimate of the true 68% and 
95% contours for the CMB. While the Ly-o constraint as used 
by WMAP pulls towards negative dn/dlnk, our Ly-a constraints 
have a strong degeneracy, with the best fit preferring less nega- 
tive dn/dlnk. Filled triangles show the models obtained from the 
CMB-only chain with dn/dlnk > 0. 

best for models with and without a set of parameters 
one wishes to test the data against. In the present case this 
reveals the success of the simplest power law scale invariant 
model with no tensors {T/S = 0, n = 1 and dn/dlnk = 0). 
The best fitted model has ~ 1449 for the combined data, 
while the best fitted model without restrictions on tensors, 
slope or running in our chains has x^ = 1445 (for the 
Ly-a forest the x^ is around 6-7 for 9 effective degrees of 
freedom, indicating a good fit to the data). The difference 
is thus only A^^ ~ 4 for effectively 3 degrees of freedom, 
corresponding to close to 68% confidence level. Adding 
these 3 degrees of freedom thus hardly improves the fit at 
all! If we fix T/S — we find that the no running model 
is within the 68% contours, with about 1/3 of all chain 
elements having dn/dlnk > 0. We also find n — 0.98 ±0.03. 
This should be compared to the WMAP results where both 
n = 1 and dn/dlnk = were exclud ed at the (roughly) 
2-sigma level in the absence of tensors llSpergel et al.ll2003l) 
(see table 8 of that paper; note that n — 1 was allowed if 
tensors were included, and that the correlation between the 
errors on n and dn/dlnk, which affects the precise level 
of confidence for the combined result, was not given). Our 
re-analysis thus reduces the statistical significance of the 



WMAP result on these two parameters. While we did not 
perform the combined analysis with 2dF it seems likely 
that with the expanded errors the addition of Ly-a does 
not improve the constraints, except for eliminating very 
strong runnings with dn/dlnk < —0.05. 

To summarize, we use the simulations described in Mc- 
Donald et al. (2003) to reanalyze the current Ly-a forest 
observations. Our results are based on hundreds of simula- 
tions that cover the range of parameter space needed with 
sufficient dynamic range. We compare simulations directly 
to the data rather than performing an inversion from the 
flux power spectrum to the linear power spectrum. We also 
emphasize the importance of external constraints such as 
the mean flux decrement, which is still poorly determined 
and can lead to a significant expansion of the errors on the 
cosmological constraints. The degeneracy direction due to 
the uncertainty in the mean fiux is similar to that from the 
CMB data. Above all, we argue that previous analyses used 
an error bar on mean flux decrement that is too small, and 
probably a value that is too low, and that better estimates 
of it and its uncertainty, combined with a grid of simulations 
that fully covers parameter space, bring the Ly-a forest data 
into better agreement with a scale-invariant model. We find 
that at the moment the combined CMB-f-Ly-a forest data do 
not rule out the simplest n = 1, T/S — and dn/dlnk = 
model. 
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